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SUMMARY 

The  non-linear  equations  of  motion  of  a spinning  projectile 
are  linearised  by  a method  different  from  that  traditionally 
used.  This  results  in  extra  terms  being  retained,  giving  a 
more  complete  description  of  the  shell  behaviour  and  revealing 
some  new  information  on  stability.  The  resulting  equations 
are  solved  as  an  eigenvalue  problem,  since  this  method  of 
solution  is  the  least  tortuous  and  produces  a good  insight 
into  the  dynamic  behaviour  of  the  shell.  It  is  found  that  the 
stability  of  the  precessional  and  nutational  modes  is  influenced 
by  gravity  and  applied  side  forces.  Results  of  numerical 
simulations  which  confirm  these  conclusions  are  given. 
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1.  INTRODUCTION 


In  recent  digital  computer  simulations  of  the  flight  of  terminally  guided 
spinning  projectiles  it  was  observed  that  lateral  control  forces  sometimes  lead 
to  processional  instability.  Since  this  result  is  not  predicted  by  the 
traditional  linear  theory  of  shell  stability,  and  to  verify  that  it  was  not  simply 
caused  by  an  incorrect  numerical  procedure,  we  have  re-examined  the  shell 
stability  theory.  It  turns  out  the  result  is  predictable  by  a linear  theory 
provided  that  the  exact  equations  of  motion  are  expressed  in  the  nonspinning  axes 
system  (p  = 0)  before  linearisation. 

A 

Nicholaides(rcf . 1)  and  Murphy(ref . 2)  have  contributed  significantly  to  the 
study  of  the  stability  of  spinning  projectiles,  and  their  results  are  still  those 
most  frequently  quoted  today.  Both  these  workers  linearised  the  equations  of 
motion,  but  they  also  assumed  that  the  flight  path  was  almost  horizontal  in  order 
to  further  simplify  the  equations.  This  meant  that  the  effects  of  gravity  and 
applied  side  forces  produced,  for  example,  by  canards  (Regan  and  Smith, 
reference  3)  did  not  enter  into  their  criteria  for  stability.  It  is  because  we 
are  here  concerned  with  quite  large  side  forces,  comparable  with  the  projectile's 
weight,  that  the  usual  approximations  become  inadequate. 

By  expressing  the  equations  of  motion  in  nonspinning  axes  and  then  linearising, 
but  without  assuming  that  the  flight  path  elevation  angle  is  small,  we  obtain  a 
system  of  six  first  order  ordinary  differential  equations,  the  solution  of  which 
reveals  the  influence  of  gravity  and  side  forces  on  stability. 

The  equations  are  solved  as  an  eigenvalue  problem.  This  elegant  technique, 
which  does  not  appear  to  have  been  used  in  shell  stability  problems  previously, 
yields  expressions  for  nutational  and  processional  damping  rates  and  frequencies 
(these  being  the  eigenvalues)  very  readily.  The  linearisation  procedure  gives 
two  more  equations  than  the  usual  treatment  does  (relating  to  the  Euler  roll 
angle  ifi,  and  the  elevation  d,  of  the  non-spinning  axes)  which  introduces  another 
mode  of  motion  which  we  shall  call  the  d mode.  It  will  be  shown  that  this 
mode  represents  the  curvature  of  the  trajectory. 

In  the  next  section  we  present  the  equations  of  motion,  first  in  the  so-called 
aerobal listic  axes  and  then  transformed  to  nonspinning  axes.  We  have  not  used 
distance  along  the  trajectory  as  the  independent  variable,  as  is  often  done. 

This  means  that  the  effect  of  drag  on  stability  is  neglected;  but  its  effect  is 
small,  and  retaining  time  as  the  independent  variable  simplifies  the  present 
analysis.  A wish  to  refrain  from  needlessly  complicating  the  treatment  with 
terms  which  produce  negligible  effect  on  stability  has  also  resulted  in  our 
excluding  Magnus  force  (however,  Magnus  moment  cannot  be  neglected).  The  effects 
of  the  neglected  terms  are  adequately  discussed  in  reference  1. 

Section  3 solves  the  linearised  equations  of  motion  as  an  eigenvalue  problem 
and  Section  4 evaluates  the  stability  criteria  from  the  derived  eigenvalues. 

It  is  found  that  the  motion  of  a typical  105  mm  shell  remains  stable  provided 
the  horizontal  component  of  the  control  force,  when  applied  at  the  nose,  remains 
less  than  about  50  N.  For  l.irger  side  forces  either  nutational  or  processional 
instability  may  be  produced  depending  upon  the  trajectory  elevation  angle  and 
whether  the  force  is  towards  the  right  or  left  of  the  trajectory.  On  the  down- 
leg,  for  example,  assuming  that  the  shell  spin  is  clockwise  when  viewed  from  the 
rear,  a horizontal  force  to  the  right,  applied  at  the  nose,  tends  to  stabilise 
precession  and  destabilise  nutation,  while  a force  to  the  left  tends  to 
destabilise  precession  and  stabilise  nutation.  The  opposite  occurs  on  the  up- 
leg  or  if  the  sense  of  shell  rotation  is  reversed,  and  the  effect  diminishes  as 
the  trajectory  elevation  angle  approaches  zero.  Instability  with  vertical 
control  forces,  applied  at  the  nose,  does  not  occur  until  the  force  exceeds  about 
1000  N.  The  effect  of  gravity  on  stability  is  very  small. 

Finally,  Section  5 gives  results  of  numerical  solutions  of  the  equations  of 
motion.  These  agree  fairly  well  with  the  analytic  solution. 


WRI:-TW-  1906(W) 


- 2 - 


2.  EQUATIONS  OF  MOTION 

The  axes  systems  used  in  formulating  the  problem  arc  illustrated  in  figure  1. 
In  addition  to  the  standard  aerodynamic  forces  and  gravity,  the  spinning  project- 
ile is  also  subjected  to  a lateral  control  force,  applied  at  the  nose,  normal  to 
its  axis,  with  components  E and  F as  indicated  in  figure  1.  Such  a force  could 
be  applied,  for  example,  by  canards  mounted  on  the  nose,  the  nose  being  mechanic- 
ally decoupled  in  roll  from  the  projectile  and  roll  stabilised  relative  to  earth. 
The  equations  of  motion  will  first  be  written  in  aeroballistic  axes,  the  angular 
velocity  of  which  is  (p  q,  r)  with  p = -r  tan  p is  the  roll  rate 

necessary  to  keep  the  y axis  horizontal.  The  x axis  lies  along  the  projectile 
spin  axis,  y lies  to  the  right  in  a horizontal  plane,  and  z in  a vertical  plane. 
These  axes  are  often  preferred  because  they  represent  quantities  in  a way  which 
is  easily  understood  by  a ground  based  observer. 

Neglecting  changes  in  spin  rate  and  axial  velocity,  which  are  small  and 
inconsequential  to  our  problem,  assuming  small  incidence,  and  putting  u = v,  the 


equations  of  motion  are  (see,  for  example,  Kolk,  reference  4): 

V = Z^v  - Vr  - r tan  O.w  + E/m  (1) 

w 3 Z^w  -f  Vq  + r tan  O.v  + g cos  0 * F/m  (2) 

q = ^q*^  “ (pA/B)  r - r tan  6 .t  - XF/B  (3) 

f = -M  V + N w + (pA/B)  q + M r ♦ r tan  0.q  XE/B  (4) 

w pw  q 


where  the  aerodynamic  force  ar.d  moment  derivatives  are  given  by: 

' M^  = p Sd.C^/BV 

^w  * »sPV^S.C^/mV 

P Sd.(qd/2V).C  /Bq 

q mq 

V “ *4  P Sd.(pd/2V).C^p^/BV 

Because  of  the  factors  (r  tan  0)  and  cos  6,  equations  (1)  to  (4)  are  non-linear. 
The  terms  involving  r tan  0 are  small,  and  a linearisation  which  neglects  them  is 
usually  a good  approximation,  particularly  if  the  flight  path  elevation  B is  small. 
However,  by  making  this  approximation,  one  obtains  a solution  whose  stability  does 
not  depend  on  the  applied  forces  and  gravity,  contrary  to  our  numerical  experiences. 
Neglecting  the  terms  containing  (r  tan  B)  is  equivalent  to  neglecting  the  residual 
roll  rate  (p  = -r  tan  B)  which  is  necessary  to  keep  the  y axis  of  the  aero- 

ballistic  axes  horizontal  as  the  projectile  pitches  and  yaws.  These  linearised 
equations  therefore  in  fact  refer  strictly  to  the  nonspinning  axes  (p  = o),  but 

A 

they  incorrectly  imply  that  E,  F and  g act  in  a constant  direction  in  this  axes 
system.  We  shall  now  show  that  if  one  starts  with  the  exact  equations  of 
motion  in  nonspinning  axes  and  then  linearises  them,  one  obtains  a set  of  equations 
which,  to  first  order,  correctly  allow  for  the  variation  in  the  direction  of  E,  F 
and  g in  these  axes. 


1 
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Transforming  equations  (1)  to  (4)  into  the  nonspinning  axes  system  is 
effected  by  rotation  about  the  x axis  by  the  angle 

where  i^l  = JT  (r  cos  + q sin  tan  6 dt)l  is  the  change  of  Euler  roll  angle 


suffered  by  the  nonspinning  axes. 

The  transformation  gives: 

V =*  2^v  - Vr  + (g  cos  + F/m)  'f>  + E/m  (5) 

w = ♦ Vq  - (E/m)  V’  - (g  sin  9^)  ® + (g  cos  6^  + F/m)  (6) 

4 = N V + M w + M^q  - (pA/B)  r + (XE/B)  ^ - XF/B  (7) 

pw  w ^ 

t * ♦ Np^^w  + (pA/B)  q + M^r  ♦ (XF/B)  ^ + XE/B.  (8) 


To  complete  this  set  of  equations  we  need  a pair  of  equations  giving  the  rates 
of  change  of  6 and  y’.  These  are  given  to  first  order  by: 

0 = q (9) 

i>  = T tan  6 (10) 

o 

Note  that  the  variables  are  now  with  respect  to  the  nonspinning  axes.  In 
deriving  equations  (5)  to  (10)  we  have  assumed  'P  is  small,  so  have  put  cos  v?  = 1, 
sin  p = p,  and  have  neglected  products  such  as  iV>  and  qy>.  The  validity  of  this 
assumption  will  be  discussed  later.  Note  that  there  is  no  loss  of  generality 
in  assuming  that  p = o initially.  We  have  also  assumed  that  the  interval  being 
considered  is  sufficiently  short  that  variations  in  6 are  small,  thereby 
replacing  6 by  (B^  + 0),  where  6 is  now  the  departure  from  its  initial  value, 

6 ; this  restricts  the  time  interval  to  a few  seconds, 
o 

It  should  be  observed  in  the  above  that  the  linearisation  of  the  transformed 
equations  is  not  equivalent  to  the  transformation  of  the  linearised  equations; 
ignoring  the  r tan  0 terms  in  equations  (1)  to  (4)  has  no  equivalent  in 
equations  (5)  to  (10).  Although  the  r tan  0 terms  are  second  order  in  magnitude, 
it  will  be  seen  that  their  cumulative  influence  upon  the  stability  of  the 
solutions  is  sometimes  comparable  to  that  of  the  first  order  pitch  damping  terms. 
This  is  a rather  sobering  example  of  the  fact  that  when  nonlinear  equations  are 
simplified,  the  relative  numerical  size  of  a term  is  not  the  sole  criterion 
determining  its  importance. 


3.  SOLUTION  BY  EIGENVECTORS 

Equations  (5)  to  (10)  give  6 linear  first  order  equations  in  six  unknowns  v, 
w,  q,  r,  p and  6.  These  can  be  written  in  matrix  form  as 

X = M X + N. 

To  determine  the  stability  of  the  solution  of  this  set  of  equations,  we  examine 
the  general  solution  of  the  homogeneous  equation 

X = M X. 


I 


A 
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Following  Struble(ref .5)  we  consider  the  linear  transformation,  u = U x,  where 

are  the  eigenvectors  of  the  matrix  M.  The  eigenvectors  satisfy  the  six 
independent  equations 

Ui  ■■  u.  , i = 1 , 6 

with  the  solutions  u^  = u^^  exp(^^t),  where  are  the  eigenvalues  of  the 
matrix  M.  They  satisfy  the  equation 

I M - t/*  11  = 0, 

where  I is  the  diagonal  identity  matrix. 

Using  the  values  for  the  elements  of  M given  by  equations  (5)  to  (10),  and  i 

expanding  the  determinant  gives  a sixth  degree  polynomial  equation  for  j 

1 

^ * a-i  ai  \l/-'  * Si,  * ai  4/  * dt  = 0 (11)  | 

where  | 

a,  - -2(mq  ♦ 2^) 

aj  = (pA/B)^  ♦ 2(M^2^  - M^V)  ♦ (M^  ^ ^o^® 

^ ^ Np(g  cos  ♦ F/m)  - (N^^E  - FM^)Z^X/Bj  g sin  O^.tan  6^. 

Tlie  other  coefficients  all  have  terms  involving  E,  F,  g,  XE  and  XF,  so  it  is  i 

simplest  to  present  them  in  tabular  form,  which  is  done  in  Table  1.  For  each  I 

entry  in  the  table,  the  dominant  component  is  given  first. 

As  an  aside,  it  is  interesting  to  note  that  for  no  external  forces  the  eigen- 
value equation  factorises  into: 

4^  - IZy,  + Mq  ♦ ipA/Blvt/  ♦ (Z^M^  - VI^)  + i|  Z^^pA/B  - VNp^^]  = 0 (12) 

and  a similar  equation  with  i replaced  by  -i,  i.e.  its  complex  conjugate;  only 
the  roots  of  equation  (12)  are  physically  admissible.  The  two  equations  give 

roots  which  are  complex  conjugates  of  each  other.  Note  that  equation  (12)  is  ^ 

basically  the  same  as  the  equation  Murphy  uses  in  his  discussion  of  projectile 
stabi lity (ref .2) . 

Returning  to  equation  (11),  the  coefficients  of  the  terms  in  the  eigenvalue 
equation  are  quite  complicated.  Fortunately  for  a tyj^ical  spinning  projectile 
a few  of  the  terms  dominate  and  the  others  may  be  neglected. 

In  Table  2 we  have  listed  representative  values  for  all  the  parameters  in  our 
problem.  Using  these  values  to  determine  the  dominant  terms,  we  obtain  the 
following  simplified  expressions  for  the  coefficients  of  equation  (11): 


a,  . -2(M^  > Z^)  (13) 
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a,  = -2(pA/B)(Z^pA/B  - VNp^)  ♦ 2VM^(M^^  + Z^)  - X(pA/B^)E  tan  0^ 


.M  = 


(15) 

(16) 


a,  = -ZVI^.g  sin  0^  - (g  sin  0^  + VZ^^)  (I^X/B)F  tan  0^ 


(17) 


a«  = M^.g^  sin  ^ 0^  * (M^  ♦ Z^^Xm/B)  (M^/m)g  sin  0^.F  tan  0^ 


(18) 


Although  interest  is  centered  on  the  eigenvalues,  for  it  is  they  that  determine 
stability,  we  will  briefly  discuss  the  eigenvectors.  They  are  given  by  the 
solutions  of  the  equations  from  which  the  eigenvalue  determinant  is  extracted. 
Since  there  are  four  unknowns  and  four  equations,  it  is  necessary  to  express 
three  of  the  variables  as  a function  of  the  fourth.  In  order  to  compare  the 
magnitudes  of  the  eigenvectors,  the  eigenvectors  need  to  be  in  the  same  units; 
so  we  will  introduce  the  dimensionless  variables  (all  of  them  angles): 


“e  ■ ^ ■ '’e/'' 


"e  • 


where  the  subscript  e denotes  that  these  are  the  eigenvectors, 
Solving  the  equations  in  terms  of  these  variables  gives; 


/J  = -i  a 
e e 


0^  - il  Z„A  - II 


S,  - I :„/»>  - II  % 


In  order  to  see  what  these  eigenvectors  look  like,  we  will  anticipate  the 
solution  of  the  eigenvalues.  It  will  be  shown  that  they  satisfy  the  inequalities 
w>  1 >X,  which  leads  to  the  approximations: 


/3. 


■1 


(-1  ♦ iZ  /cJ)  a 
V yi'  ' e 


Typical  values  of  Z^/co  are  0.03  for  precession  and  0.0013  for  nutation,  from 

which  it  follows  that  the  angular  excursions  of  the  eigenvectors  are  of  nearly  the 
same  magntiude,  and  they  are  close  to  phase  quadrature.  When  the  motion  is 
constrained  to  rotation  about  the  center  of  gravity,  i.e.,  Z^  = 0,  the  eigen- 
vectors are  exactly  in  phase  quadrature.  It  is  seen  that  having  no  such 
constraint  makes  only  a small  difference  to  the  eigenvectors. 

Addition  of  gravity  and  side  forces  complicates  the  eigenvectors  enormously. 

Not  only  are  the  expressions  for  modified,  but  there  are,  of  course, 

two  extra  terms,  corresponding  to  ip  and  0.  However,  the  simple  relationships 

I 
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6 = TJ 

e c 


'P  = tan  0 f 
e o e 


show  that  the  >fi  eigenvector  vanishes  for  horizontal  flight  but  the  0 eigenvector 
does  not. 


4.  STABILITY  CRITERIA 

The  frequency  and  damping  of  the  normal  modes  of  the  projectile's  motion  are 
given  by  the  imaginary  and  real  parts  of  the  eigenvalues  of  equation  (11). 

Rather  than  attempt  to  attack  this  equation  directly,  we  shall  adopt  an  approach 
which  makes  use  of  the  known  properties  of  the  solutions. 

As  discussed  earlier,  when  gravity  and  side  forces  are  not  present,  the  eigen- 
value equation  reduces  to  a fourth  degree  one  whose  solution  is  well  known, 
consisting  of  two  pairs  of  complex  conjugates.  The  motions  corresponding  to 
these  are  the  precessional  and  nutational  modes; 

'^P  = 

}j/  = X ± iw 

N N N 

whose  relative  magnitudes  for  a typical  spinning  projectile  are: 

C^>U^>1,  ixi^  'l\lp<1. 

Introduction  of  equations  (9)  and  (10)  increases  the  degree  of  the  eigenvalue 
equations  by  two  and  so  introduces  two  more  roots,  which  we  shall  denote  by 
^0  These  roots  may  either  be  a complex  conjugate  pair,  or  both  real.  For 

now  we  shall  assume  they  are  of  the  form 

'^(9  ,<p  = ^0  H 

where  ^ is  imaginary  when  the  roots  ^ are  real.  The  roots  tend  to  zero 

for  vanishingly  small  values  of  g sin  0,  E and  F,  and,  as  will  be  shown,  they  are 
always  much  less  than  unity. 

Expressing  the  roots  in  this  form  gives  a sixth  degree  polynomial  for  the 
eigenvalues : 

^ ♦ bi  ^ bj  ^ + bs  ^ + b4  + bj  ^ + b6  =0  (19) 

where  the  coefficients  are  sums  and  products  of  the  X's  and  cJs.  By  retaining 
only  the  dominant  terms  these  simplify  to 

b.  = (20) 

b,  = 


(21) 
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(22) 

b,  . b,,> 

(23) 

(24) 

(25) 

We  ar'  now  in  a position  to  obtain  explicit  expressions  for  the  eigenvalues  in 
terms  of  the  parameters  of  the  problem.  First  we  confirm  that  ^ and  ^ are 

indeed  small;  from  equations  (16),  (17),  (18),  (23),  (24)  and  (25)  we  have 

= 8 ^ ^^w  ^ 8 sin  6JV)  (X/2VM^B)  F tan  6^  (26) 

= I (M^^  - Xm  g sin  d^/BV)/(V^  M^^tn)!  8 sin  tan  9^  (27) 

The  largest  term  in  either  of  these  expressions  is  g sin  0 /V,  which  is  of  the 
order  of  0.04. 

Equations  (14)  and  (21)  give  immediately  for 

= pA/B  - VM^B/pA  (28) 


correct  to  first  order  in  the  binomial  expansion  for  the  square  root.  We  can 
now  calculate  from  equations  (16)  and  (23): 


Wp  = VM^B/pA  + M^(B/pA)^  (29) 

Also,  from  equations  (13),  (15),  (20),  (22)  we  have: 

X^  = H 1 + 2VM^^(B/pA)M  VNp^  B/pA 

(Mq  - Z^^)  (B/pA)'  VM^  - (X/2pA)  |l  + 2VM^(B/pA)ME  tan  0^.  (30) 

Finally,  from  equations  (15)  and  (22): 

Xp  = - (1  + 2VM^^(B/pA)M  VNp^B/pA 

- (M^  - Z^)  (B/pA)'  VM^  + (X/2pA)  | 1 + 2VM^(B/pA)ME  tan  9^ 

- g sin  9 /V  - (Z  + g sin  0 /V) (X/2VM^B)F  tan  9 . (31) 

o w o w o 


! 


4 


In  every  case  we  have  put  the  dominant  terms  first  in  the  expressions  on  the 
right  hand  side.  When  E = F = g = 0 these  expressions  reduce  to  the  well  known 
solutions  of  the  linearised  equations  of  motion(ref . 1) . Pitch  damping  and  normal 

force  are  the  prime  stabilising  factors  for  nutation  and  precession  respectively. 
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If  the  Magnus  moment  is  large,  one  or  other  of  precession  and  nutation  becomes 
unstable.  The  precessional  and  nutationaJ  frequencies  are  seen  to  be  independent 
of  gravity  and  side  forces  to  second  order,  whereas  the  damping  depends  on  these 
forces  to  the  first  order  of  approximation.  However,  the  terms  involving  the 
vertical  applied  force,  F,  arc  generally  much  smaller  than  those  involving  the 
horizontal  force,  E.  For  the  time  being,  therefore,  we  shall  neglect  the  terms 
in  F . 


A given  mode  is  unstable  if  the  corresponding  A is  positive.  From  the  sign 
of  each  term  in  the  expressions  for  A we  can  ascertain  whether  given  forces  are 
stabilising  or  destabilising.  This  information  is  summarised  in  Table  3.  For 
both  nutation  and  precession,  the  destabilising  force  must  be  greater  than  a 
certain  minimum  value  to  overcome  the  stabilising  influence  of  M and 


respectively. 

Substituting  values  for  the  parameters  listed  in  Table  2 gives  the  following 
approximate  expressions  for  the  eigenvalues: 


w = 4.5  rad  s ' 

P 

A = -0.269  + 0.0067  E tan  0 - 0.039  sin  0 s 

p o o 

= 105  rad  s’* 


A„  = -0.353  - 0.0067  F.  tan  0 s 

N o 


The  contribution  of  the  vertical  force,  F,  has  again  been  omitted  as  it  is 
smaller  than  the  other  terms. 

We  have  plotted  damping  rates.  A,  in  figure  2.  It  is  seen  that  the  projectile 
is  nutationally  unstable  for  E tan  0 < - nON,  and  precessionally  unstable  for  E 
tan  0 > 40N. 

The  9,  ifi  mode  is  a mode  of  motion  which  arises  from  the  introduction  of 
equations  (9)  and  (10),  which  allows  0 and  'fi  to  vary.  Note  that  this  mode  has 
a time  constant  of  the  order  to  25  s and  so  does  not  vary  significantly  during  the 
few  seconds  of  validity  of  the  solution.  Its  physical  significance  is  that  it 
represents  the  curvature  of  the  trajectory.  We  shall  demonstrate  this  for  the 
case  where  gravity  is  the  only  external  force. 

First,  let  c be  the  elevation  of  the  tangent  to  the  trajectory.  By  relating 
its  rate  of  change  to  the  force  normal  to  the  trajectory  we  have 


de/dt  = -g  cos  c/V. 


Now  e =9-0,  where  a is  the  vertical  component  of  incidence,  positive  when  the 
shell  nose  is  above  the  tangent  to  the  flight  path.  Replacing  9 by  9^  + 9, 

expanding  cos(9  - a),  and  using  the  small  .ingle  approximations  sin  o = a, 
sin  9 = 9 , gives 

d(9  - a)/dt  = (g  sin  9^/V)  (9  - n)  - g cos  9^/V 

The  general  solution  of  this  equation  has  a time  constant  V/g  sin  9^,  which  is 

the  same  as  for  the  eigenvalue  for  the  9,  yj  mode. 

Second,  since,  from  equation  (5), 

r ~ (g  cos  'f>  * .. . 


1 

! 

i 

i; 

I 


i 

i 

I 


ijM 


9 


WRF.-TR-1906(W) 


i 


i 


i 


it  follows  from  equation  (10)  that 

^ ~ (g  sin  ^q/V)  * ...  , 


the  solution  of  which  also  has  the  time  constant  V/r  sin  0 . 

o 

This  confirms  our  interpretation  of  the  6,  <p  mode.  The  degeneracy  of  the 
roots  of  the  6,  <fi  mode  is  removed  by  the  normal  force,  F.  The  form  of  the  roots 
of  the  6,  tp  mode  depends  on  tho  sign  of  ^ given  by  equation  (27).  Since 

I - Xmg  sin  0^/BVl  is  positive  for  the  parameter  values  given  in  Table  2,  the 

sign  ^ ^ ^ then  the  same  as  the  sign  of  F.  For  F positive,  i.e.  downward 

force,  the  two  eigenvalues  are  a pair  of  complex  conjugate  niunbers.  For  F 
negative,  we  put  ^ a eigenvalues  have  the  real  values, 

X ± AX. 


5.  NUMERICAL  COMPUTATIONS 

There  is  always  satisfaction  when  exact  numerical  calculations  confirm  an 
analytic  treatment.  This,  plus  the  benefit  of  computer  plotted  figures  showing 
the  nutational  and  precessional  instabilities  very  clearly,  arc  our  excuses  for 
presenting  the  results  of  some  numerical  simulations  in  this  section. 

We  have  written  a 6 degree  of  freedom  computer  program  to  calculate  the 
projectile  motion  numerically.  The  equations  solved  are  (1)  to  (4)  plus  the 
equations  of  motion  in  the  x direction: 

u * S P ■ 8 ^ 

i>  - S P V^S.C^p.pd/2V 


with  C = -0.13  and  C,  = -0.025.  We  have  also  added  a Magnus  force  term  to 

X X p 


equations  (1)  and  (2),  with  C^  = -1.5. 

TpU 


By  running  the  program  with  and 


without  these  refinements,  the  contribution  of  these  added  refinements  was  found 
to  be  small,  but  they  have  been  added  for  completeness. 

The  results  of  three  runs  of  the  program  are  given  in  figures  3 to  5.  The 
shell  parameters  were  set  at  the  values  given  in  Table  2,  which  apply  to  halfway 
down  the  downleg  of  a typical  105  mm  shell  trajectory.  The  elevation  was  set 
initially  at  -43^  and  the  velocity  at  250  ms  ' . Ouring  the  4 s for  which  the 
trajectory  was  calculated  the  mean  elevation  decreased  to  -48°  and  the  velocity 
increasea  to  264  ms  ' ; i.e.  these  values  did  not  depart  too  much  from  the 
constant  value  assumed  in  the  analytic  solution.  Tlie  numerical  calculations  used 
a standard  fourth  order  Runge-Kutta  routine.  The  time  step  was  0.0025  s;  if  it 
was  much  larger  than  this  nutation  was  found  to  be  artificially  damped. 

Each  figure  shows  the  incidence  (0,  a)  and  the  angular  rate  (q,  r) ; it  also 
shows  values  of  "P  and  its  integral  <p,  the  Euler  roll  angle  of  the  nonspinning 
axes,  as  functions  of  time.  These  figures  show  very  clearly  the  conditions  for 
instability  which  were  discovered  analytically.  With  no  side  forces  acting 
(figure  3),  both  precessional  and  nutational  modes  are  stable;  a 100  N force  to 
the  right  (figure  4)  produces  nutational  instability  and  increases  the  precessional 
stability,  whereas  a 50  N force  to  the  left  (figure  5)  gives  precessional 
instability  while  increasing  nutational  stability. 
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Tho  processional  and  nutational  modes  are  easily  distinguished  in  figures  3 to 
5,  so  estimates  can  be  made  of  their  frequency  and  damping.  The  results  are 
also  plotted  in  figure  2.  Comparison  with  the  eigenvalues  of  the  analytic 
treatment  shows  a systematic  difference  in  X of  about  0.01  s due  to  the  neglect 
of  drag  in  tho  analytic  treatment.  Murphy(ref . 2)  has  shown  that  drag  has  a 
destabilising  effect,  increasing  X by  'j  p V^SC^/mV.  For  the  case  we  are 

considering  here,  this  equals  0.01  s~‘ . In  addition  to  the  systematic  difference, 
there  is  a scatter  of  about  0.005  s"' , which  is  largely  due  to  errors  in 
estimating  X from  the  numerical  simulations. 

Finally,  the  numerical  calculations  enable  us  to  see  how  valid  was  our 
assumption  that  the  integrated  residual  roll,  'P , was  small.  It  is  seen  that 
although  p is  very  small  0.01  rad)  when  there  are  no  side  forces,  the 
approximation  is  only  valid  on  short  time  intervals  for  large  side  forces,  where 
p can  approach  1 rad  after  a few  seconds. 


6 . SUMMARY 

The  equations  of  motion  of  a spinning  projectile  have  been  written  in  non- 
spinning axes,  whose  angular  velocity  is  (0,  q,  r) , in  which  the  directions  of 
gravity  and  side  forces  are  not  constant,  but  wobble  in  the  complementary  manner 
to  the  Fulerian  roll  angle  of  the  nonspinning  axes  system.  By  assuming  that 
this  wobble  is  small,  an  analytic  solution  of  the  linearised  equations  is  obtained 
whose  general  solution  contains  gravity  and  side  forces  in  the  damping  term. 

It  is  found  that  on  the  trajectory  downleg,  a horizontal  side  force  to  the 
right,  applied  at  the  nose,  tends  to  stabilise  precession  and  destabilise  nutation, 
and  a ‘‘orce  to  the  left  tends  to  destabilise  precession  and  stabilise  nutation. 

The  opposite  occurs  on  the  upleg.  This  means  that  a large  enough  side  force 
(greater  than  about  40  N,  for  a 105  mm  shell,  when  applied  at  the  nose)  will 
always  cause  the  motion  of  the  projectile  to  be  unstable. 

Cravity,  and  side  forces  in  the  vertical  plane,  have  a much  smaller  effect  on 
precessional  and  nutational  stability  than  do  horizontal  side  forces.  For 
expected  parameter  values,  these  forces  are  unlikely  to  be  large  enough  to  over- 
come  the  stabilising  effect  of  the  aerodynamic  damping  moment  and  normal  force. 

Nunerical  simulation  of  the  motion  of  a projectile  using  a 6 degree  of  freedom 
computer  program  has  confirmed  these  predictions  of  the  analysis. 
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NOTATION 

coefficients  in  eigenvalue  polynomial  equation 
moment  of  inertia  about  x axis 

coefficients  in  eigenvalue  equation  in  terms  of  roots 

moment  of  inertia  about  y and  z axes 

roll  damping  moment  coefficient  derivative 

pitching  moment  coefficient  derivative 

pitch  damping  moment  coefficient  derivative 

normal  force  coefficient  derivative 

Magnus  moment  coefficient  derivative 

axial  force  coefficient 

Magnus  force  coefficient 

reference  length  (=  projectile  diameter) 

applied  side  force  in  horizontal  direction 

applied  side  force  in  vertical  direction 

gravitational  constant 

identity  matrix 

mass 

matrices  in  Section  3 

pitching  moment  derivative  (equation  (3)) 
pitch  damping  moment  derivative  (equation  (3)) 

Magnus  moment  derivative  (equation  (3)) 
projectile  roll  rate 

roll  rate  of  axes  system  ( = -r  tan  B in  aeroballistic  axes 

L = 0 in  nonspinning  axes 

angular  velocity  about  y axis 

angular  velocity  about  z axis 

reference  area  (=  »d^/4) 

time 

eigenvector  of  motion 

velocity  component  in  x direction 
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V 

V 
w 

X 

X 

y 

z 

X 

z 

w 

Greek  letters 
a 

e 

X 

0 

e 

o 

p 

<jj 


aeroballistic 

axes 


velocity  component  in  y direction 
projectile  velocity  (in  still  air) 
velocity  component  in  z direction 
vector  (v,  w,  q,  r,  >P,  0) 
along  projectile  ^ 
horizontal 
in  vertical  plane 
distance  of  side  force  forward  from  c.g, 
normal  force  derivative  (equation  (1)) 

vertical  component  of  incidence,  positive  nose  up 

elevation  of  tangent  to  trajectory 

damping  (real  part  of  eigenvalue) 

elevation  of  the  x axis 

initial  elevation 

atmospheric  density 

residual  roll  angle  of  nonspinning  axes  (=  -/  r tan  0 dt) 

eigenvalues  of  motion 

frequency  (imaginary  part  of  eigenvalue) 


Subscripts  and  superscripts 


N 

P 

0,  ^ 
a 

a 


nutational  mode 
precessional  mode 
, 'P  mode 

differentiation  with  respect  to  time 
matrix  or  vector 
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TABLE  2.  TYPICAL  VALUES  FOR  PARAMETERS 


Primary  parameters 

Derived  parameters 

V 

250  ra.s”' 

P 

1050  rad.s  ‘ 

P 

1.05  kg.m  ^ 

m 

15  kg 

A 

0.023  kg.m^ 

B 

0.22  kg.m^ 

d 

0.105  m 

S 

0.0087  m^ 

X 

0.3  m 

3.8 

M 

2.07  m"'.  s ' 

Sia 

-1.8 

Z 

w 

-0.137  s'* 

C 

-16 

M 

-0.457  s'' 

mq 

q 

0.048  m"'.  s'' 

c 

0.4* 

N 

npa 

pw 

E 

0 * 100  n"^ 

F 

0 •*  100  n"^ 

g 

9.8m.  s"* 

The  values  apply  to  halfway  on  the  downleg  of  a 105  mm  shell  trajectory. 

* Magnus  moment  increases  rapidly  from  0.02  to  1.0  in  the  transonic  region. 

We  have  taken  a mean  value. 

+ Canards  mounted  on  the  nose  are  likely  to  produce  side  forces  in  this  range. 


TABLE  3.  EFFECT  OF  FORCES  ON  NORMAL  MODES 


Force 

Precession 

Nutation 

g sin  6 > 0 

stabilising 

nil 

< 0 

destabilising 

nil 

E tan  6 > 0 

destabilising 

stabilising 

< 0 

stabilising 

destabilising 

F tan  6 > 0 

stabilising 

nil 

< 0 

destabilising 

nil 

I 6 is  positive  on  the  upleg 

E is  positive  to  the  right  viewed  from  the  rear 
F is  positive  downwards 

Projectile  spins  clockwise  as  viewed  from  rear 
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Figure  1.  Axes  systems 


Side  force  Etan^  N 


= analytic  solution 

= from  numerical  simulation 

These  calculations  apply  to  half  way  on  the  downleg  of  a 105  mm  shell 
trajectory. 

Figure  2.  Comparison  of  analytic  solution  with  numerical  simulation 
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